Localization of charged quantum particles in a static 
random magnetic field 



A.G. AronovQ, A.D. Mirliifl and P. Woffle 

Institut fur Theorie der Kondensierten Materie, 
Universitdt Karlsruhe, 76128 Karlsruhe, Germany 
(December 20, 1993) 

Abstract 

We consider a charged quantum particle in a random magnetic field with 
Gaussian, delta-correlated statistics. We show that although the single par- 
ticle properties are peculiar, two particle quantities such as the diffusion con- 
stant can be calculated in perturbation theory. We map the problem onto a 
non-linear sigma-model for Q-matrices of unitary symmetry with renormal- 
ized diffusion coefficient for which all states are known to be localized in d = 2 

dimensions. Our results compare well with recent numerical data. 
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The problem of a charged quantum particle moving in a static random magnetic field in 
two dimensions has received renewed interest recently For one, the problem may be 

considered as a limiting case of a system of particles interacting via a gauge field. Models 
of this type have been proposed to describe a state with charge-spin separation of the 
conduction electrons in High-T c superconductors |3||| . Secondly, an experimental realization 
of a random magnetic field due to the pinned vortex lines of a superconducting layer on top 
of a semiconductor heterostructure has been reported recently f7|. Thirdly, the problem is 
thought to be relevant for the quantum Hall effect in the limit of the half-filled Landau level 



m 



A number of numerical investigations have been performed, with conflicting results. In 
Ref. it was argued on the basis of results of numerical diagonalization on square lattices 
of up to 10 4 sites for zero- average random flux per plaquette and in addition site-diagonal 
disorder and a uniform magnetic field that localization can be suppressed by the random 
flux. In Ref. [|J the conductance of a square lattice of quantum wires subject to an random 
magnetic flux per plaquette, distributed uniformly between —(fio/2 and <fio/2 was calculated 
numerically (0o is the A ux quantum). Although no definite conclusion could be drawn, the 
results were found to be consistent with the existence of extended states and a mobility 
edge. In contrast, the results of applying the finite size scaling method of MacKinnon and 
Kramer to the random magnetic field problem reported in Ref. [f|] suggested that all states 
are localized by a random magnetic field. Since the localization length for a two-dimensional 
disordered system may be very large, it is obvious that numerical studies of systems of finite 
size are of limited value in deciding the principal question whether there exist extended 
states in these systems. 

In this letter we show that the problem of charged particles in a static random magnetic 
field can be mapped onto a nonlinear sigma-model of unitary matrices. The latter model 
has been proposed for disordered systems containing random spin scattering centers as well 
as models featuring random phase fluctuations of the hopping matrix elements of a tight- 
binding Hamiltonian [0,0]. Perturbation theory for this model yields a divergent quantum 



correction to the conductance in two-loop order ||14|| . As a consequence the scaling function 
in two dimensions remains negative, leading to the result that all states are localized for 
these models. 

We study the transport properties of a charged spinless quantum particle (mass m, charge 
e) in two dimensions in a static random magnetic field H(r) = V x A(r) normal to the plane, 
as defined by the Hamiltonian 

H=^~— P-A + ^—J 2 , (1) 
2m mc 2mc z 

where p = —iW is the momentum operator and the Coulomb gauge (V ■ A = 0) has been 
used (h = 1). The magnetic field is assumed to be Gaussian distributed and delta-correlated, 
with vanishing mean and variance in Fourier space 

(H(q)H(-q)} = (h 2 ) (2) 

Accordingly, the variance of the vector potential is (A a (q)Ap(—q)) = (l/q 2 )(h 2 ) S^Jq) where 
^apio) = <W — q&% accounts for the transverse character of the field (q — qj \ q |). By 
contrast, an independent random distribution of phases of the hopping matrix elements 
correspons to a delta-correlated distribution of vector potentials given by (A a (q)Ap(— q)) = 
(a 2 )^J /3 (g). We will comment on this case later. 

Assuming the fluctuation strength of the magnetic field, (h 2 ), to be weak, we first consider 
perturbation theory. The standard Feynman diagram language for impurity scattering may 
be employed, with the impurity line describing scattering of a particle from momentum state 
| p + q/2) into state | p ' + q/2) and a hole from | p — q/2) into | p ' — q/2) being given by 

Wpp ,{q) = v 2 Q k- 2 {{p + p'f -q 2 - [{p + p') ■ k] 2 + (q- k) 2 } (3) 

where k = p — p' is the transferred momentum and v 2 = e 2 (h 2 )/4m 2 c 2 is a velocity squared 
characteristic of the strength of disorder. Note that w^* is strongly singular in the forward 
direction (p = p'), due to the long range of the vector potential fluctuations, even though the 
magnetic field fluctuations are assumed to be short-ranged. In lowest order the imaginary 



part of the (advanced) single particle self energy S (p, E) on the energy shell (E = p 2 /2m) 
is given by 



2r ■/ zm 2m 



Here (dp) = d 2 p/(2n) 2 , p ■ p' = cos (ft and iVo = f^- is the density of states. The integral 
in (H) is strongly divergent at = 0, which may be traced to the contribution of vector 
potential fluctuations in the limit q — > 0. A self-consistent treatment of the divergence 
leads to a much weaker dependence as discussed below after (|23|). Nonetheless, the single 
particle relaxation rate l/2r can be expected to diverge, which might lead to the generation 
of a branch cut in the single particle Green's function, as argued in Ref. [T||J] (see however 
|18|). We will show below that the transport relaxation rate and the diffusion constant are 
nonsingular, in agreement with results of a naive perturbation theory || . This may be done 
by regularizing the infrared divergence in (f|) in a convenient way. We will use a soft cut-off, 
replacing cot 2 | in ([|) by w{(f) = cos 2 |/(sin 2 | +7 2 ), although the precise form of the cutoff 
is not important. At the end of the calculation we take the limit 7 — > 0. 

Next, let us consider the diffusion propagator T ( "diffuson" ) , obtained by summing the 
particle-hole ladder diagrams 



rV(g» = v) ,@ +J(dp) w^($G* : (E + )G^(Ejr p/ , ff (q,u) (5) 

where 

Gf\E) = \E-^±—V (6) 

p v ; L 2m 2rJ v ; 

and p± = p" ± q/2, E± = E ± u/2. The solution of (^) may be easily obtained in terms 
of the eigenfunction of the operator w^>(q = 0) (see Ref. ]15|). In two dimensions and for 
I p I = I p ' I one has 

Wpp'(q= 0) = W nX*n(P)Xn(p ') (7) 



with 



w* = vl r %{<P)e~^. (8) 
JO Ztt 



and 



Xn(p) = e m * (9) 
Here 4> is the polar angle of p. The leading contribution to Tppi is obtained as 

V ^' M= -^Dq^N^ (10) 
where the diffusion coefficient D is given by 

D=\v\{l + -^) = l -v\ tr (ff) 
2 wo — wi 2 

and v = p/m. Here we have defined the transport relaxation time r tr as 

1 2 [ d< P ,2 0/! J.\ 2 / 10 x 
= mi) / — COt — (1 — COS 0) = 77W , (12) 

r tr J 2n 2 

In contrast to r, the transport time r tr is finite in the limit 7^0. 

In a time reversal invariant system, the coherent backscattering described by a diffusion 
pole in the particle-particle ladder diagrams, the socalled Cooperon, plays a dominant role. 
The Cooperon Cpp>(q,uj) obeys the integral equation fl5|), with wppi(q) replaced by wtti{Q), 
where k = \{p — p' + q), k ' = \{p' — p + q) and Q = p + p' . The fact that the vector 
potential A couples to the particle momentum p (see (|l|)), which changes the sign under 
time reversal, leads to the relation Wtti{Q) = —wpp '(([)■ Correspondingly, the Cooperon is 
finite in the limit q, u — > and cannot play any role in bringing about localization in the 
present case. 

In the following we map the problem onto a nonlinear sigma-model of unitary symmetry. 
The generating functional for two-particle Green's function of the retarded-advanced type 
(RA) may be represented in terms of a functional integral over a supersymmetric field 



if) = (tpi, xi, f2, X2), where <p\p are bosonic and xi,2 are fermionic components as jTI 



where 



Z = J Dty]exv-(S + Sx) (13) 



S = -ij d 2 r{i,A(E + ^V 2 )^ - (14) 



and 



S x = -z JL J d 2 r{^A{-iA ■ V)^}. (15) 

The 4x4 matrix A is diagonal, A = diag(l, 1, — 1, —1). 

Averaging over the vector potential one finds that Si in (15) has to be replaced by 

si ff = -i e Kn^r&^s^ ( 16 ) 

Q<10 

where w^^, is defined in (|J), and only long- wave length fluctuations are considered (q < qo). 
It is useful to introduce the representation of wgg, in terms of eigenfunctions (^), and to 



define "density" fields |L5| 



PZ=E^I + ^Xn(k) (17) 
k 

in terms of which Si can be written as 

si n = -*[ d 2 r y: wrfifwdrw- (is) 

As usual, the interaction term may be decoupled with the aid of Hubbard- Stratonovich fields 
Q n (r), which are (4 x 4) supermatrices of unitary symmetry. The functional integration over 
the primary fields may be performed, yielding 



£>[Q]exp-S{Q} (19) 
where the effective action of the Q-fields is at first given by 

S= />r{Str In G _1 - - £>„ Str Q 2 } (20) 



where Str denotes the supertrace and 



Gp 1 - e 



p 

2m 



+ iJ2 W nXn(P)Qnk 



(21) 



and e = diag (E + f , E + f , E - f , E - f ). The saddle point of exp(-S) is at Q n = Q$, 
where is 



Q [ n ] = * I (d P )Xn(p)G(p)A 



i5 n0 / (dp) 



<- - — + iw Q^A 
2m 



so that the Green's function at the saddle point is given by 

-2 



A 



(22) 



G(p) 



p i . 

e — — 1 A 

2m 2t 



(23) 



Equations fl22"|) and (^) are a statement of the self-consistent Born approximation for 
the single-particle relaxation rate -. In contrast to the lowest order expression (^) for 
-, which scales with the cutoff 7 as r _1 oc 7 _1 , the self-consistent value is given by 
t -1 = [^P 2 Vq In (l/£'r7)] 1 / 2 and hence shows a weaker divergence as 7 — ► 
We now expand the action around the saddle point: 



S = S + ^J {dq){ J {dp)Y, w nXn{P)w m Xm{p + q) Str G{p)A5Q n {q)G(p + q)A5Q m (-q) 



-J2w n Str 5Q n (q)5Q n (-q) } 



(24) 



The p- integral is only finite for the products G R G A , which are generated by the off-diagonal 
components of 5Q, denoted 5Q. In the limit of small q,u, using 



/ (dp) X n(p)Xm(p)G A (E - ^,P)G R (E + + q) 



= 27iN r\5 nm + ir(u5 nm 

L v m 

where 6$, = (XnXm cos £ 0), ^ = 1,2, we find 



pqr\2 



m 



(25) 



5 = S n 



> n (l--^)Str [5Qn(g)5Q„(-g)] 

W7q 



^ J n 

+w r(-iu + D q 2 ) Str [SQ {q)SQ {-q)] 
£ Str[5g (g)5Q„(-g)]} 



2 m 



(26) 



n=±l 



As expected, the n = mode is massless and describes interacting diffusons. The coeffi- 
cient of th g 2 -term is the bare diffusion constant, D = \v 2 t, which tends to zero if the 
infrared cutoff is taken to zero. However, the bare diffusion constant gets dressed by the 
coupling to the massive n — ±1 modes. Indeed, integrating out 5Q\ and SQ^i produces a 
renormalization term 



AS = ~f (dq)D r-^^q 2 Str [6Q {q)6Q Q (-q)] (27) 



which combined with the bare diffusion term has the effect of changing the bare diffusion 
constant D into the renormalized D as defined in fll2j) . The final result for the effective 
action is 

S = Sq 7T-/Vq J (dq){(Dq 2 - iu) Str [5Q (q)5Q (-q)}} (28) 

We note in passing that the coefficients of terms with higher spatial derivatives of Q will be 
renormalized in a similar way and may be expected to be finite as well. 

The expansion ( pip of the action in terms of 5Qo serves to determine the coefficients 
of the two terms in the nonlinear sigma model obtained from fl2"0|) by keeping only the 
integration over the saddle point manifold: 

S„ = ^ Jd 2 r{-D Str ( VQ • VQ) - 2iuo Str (AQ) } (29) 

where the rescaled field Q(r) is constrained by Q 2 (r) = 1 . 

We have thus shown that the problem of a charged quantum particle in a static random 
magnetic field is equivalent to a nonlinear cr-model of interacting Q-matrices with unitary 
symmetry. This model has been studied extensively []12-14|. It is known that the Gell- 



Mann-Low /3-function describing the scaling behavior of the dimensionless conductance g 
(in units of e 2 /h) with the length of the sample is given inleading order for large g b 

It follows then that all states are localized, and that the localization length in the weak 
disorder regime is given by 
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£ = £oexp(7r 2 £ 2 ) (31) 

where g = mv 2 T tr /2 = (v 2 /4vq) is the Drude conductance, and £ — VT tr- For the case of a 
delta-correlated distribution of vector potentials the single particle relaxation rate does 
not show an infrared divergence, and the nonlinear sigma model may be derived in the usual 
way. The more complete derivation given here leads to a renormalization of the diffusion 
constant as D = 2D = (47rN e 2 (a 2 ) /c 2 ) -1 where (a 2 ) was defined after (2). 

Our results are in good agreement with the available numerical data |||| . The authors of 
these papers studied the lattice version of the problem with the maximum possible disorder 
corresponding to Vq ~ v. Accordingly, the typical values of the "bare" conductance go are of 
order of unity. However, when approaching the center of the band (i.e. when g increases), 
the localization length was found to grow exponentially Q, in agreement with (31). The 
finite size scaling analysis yields for the localization length A of a quasi ID strip of width 
M: 

A(M) = Mf{t/M) (32) 

Calculating the scaling function f(x) by using (31) and comparing with the known result 
for the localization length of a quasi ID system [|16||, we obtain [|17||: 



!-\/ln x, x ^> 1 
(33) 
x, i C 1 

This agrees well with the asymptotic behavior of f(x) as obtained by numerical means in [|J, 
for both x ^> 1 and x C 1 (see Fig.l). Scaling behavior of the conductance g(L/£) obtained 



in || (see Fig. 4 of Ref. 0) is also compatible with the scaling law g(L) ~ (l/7rWln(£/L) 
(f > L) which follows from Eqs.flSIJ), fl£Tj). 

To summarize, we have shown here that the nonlinear sigma model description of disor- 
dered systems with broken time reversal invariance holds true even in the case of long-ranged 
fluctuations of the vector potential when the single particle properties are dominated by in- 
frared divergencies. This is true provided the magnetic field fluctuations are short ranged. 
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In the opposite case of long-ranged magnetic field fluctuations it is conceivable that even the 
transport relaxation rate diverges, signalling a different physical regime. One may speculate 
that then the topological excitations governing the behavior in the Quantum-Hall- Effect, 
where the average magnetic field is finite and large, will play a role. However, for short- 
ranged magnetic field fluctuations the topological term is absent. Finally, we emphasize that 
our analysis is restricted to quenched random magnetic fields. To what extent this model 
applies to dynamical gauge field models remains to be seen. 

This work was supported by the Sonderforschungsbereich 195 der Deutschen Forschungs- 
gemeinschaft (A.G.A.) and by the Alexander von Humboldt Foundation (A.D.M.). We thank 
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FIGURE CAPTION 



Fig.l. The scaling function /(x), Eq. (p2|) , as obtained by numerical study of the problem 
in Ref. M . The dashed and dot-dashed lines represent the asymptotical behavior for i < 1 
and x ^> 1 respectively, given by Eq. (|33|) . 
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